Scaling behavior of impurities in mesoscopic Luttinger liquids 
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Using a functional renormalization group we compute the flow of the renormalized impurity potential 
for a single impurity in a Luttinger liquid over the entire energy range - from the microscopic scale 
of a lattice- fermion model down to the low-energy limit. The non-perturbative method provides a 
complete real-space picture of the effective impurity potential. We confirm the universality of the 
open chain fixed point, but it turns out that very large systems (10*- 10^ sites) are required to reach 
the fixed point for realistic choices of the impurity and interaction parameters. 



The low-energy physics of one-dimensional interacting 
electron systems with Luttinger liquid (LL) behavior is 
dramatically affected by the presence of a single impurity. 
10-^ The problem is usually mapped onto an effective 
field theory using bosonization, where terms which are 
expected to be irrelevant in the low-energy limit are ne- 
glected. [0-^ Then forward and backward impurity scat- 
tering decouple, and the more important backscattering 
processes are modeled by a single amplitude Vb ■ From a 
perturbative bosonic renormalization group (RG) calcu- 
lation and a boundary conformal field theory analysis 
the following picture emerged: In a chain of spinless 
fermions j^] with repulsive interactions (LL parameter 
Kp < 1) the backscattering amplitude Vb is a relevant 
perturbation which grows as A^*""^ when the flow pa- 
rameter A is sent to zero, and the perturbative analysis 
breaks down. This behavior can be traced back to the 
power-law singularity of the 2kF density response func- 
tion in a LL. ||^,^ On the other hand a weak hopping tw 
between the open ends of two semi-infinite chains is ir- 
relevant and scales to zero as A^p ^^-0 Assuming that 
the open chain represents the only stable fixed point it 
was argued that at low energy scales and even for a weak 
impurity physical observables behave as if the system was 
split in two chains with open boundary conditions at the 
end points. [Q Here we focus mainly on the local spec- 
tral weight pj{uj) for lattice sites j close to the impurity 
and energies u! close to the chemical potential /i. For 
Pj{oj) a power-law suppression Pj{uj) ^ \lu\°''^ with the 
boundary exponent as — — 1 which only depends 
on the interaction strength and band filling, but not on 
the impurity parameters, was predicted. Within the 
bosonic field theory the above conjecture was verified by 
refermionization quantum Monte Carlo calculations, 
[^,|| and the thermodynamic Bethe ansatz. |0 

To confirm the field theoretical scenario and the va- 
lidity of the underlying assumptions for a microscopic 
fermionic system with LL behavior, numerical methods 
(exact diagonalization [ED] , density-matrix renormaliza- 
tion group [DMRG]) were applied to the lattice model 
of spinless fermions with nearest neighbor interaction. 



lH'O 1^1 Comparing ED data for up to = 23 sites 
with the field theoretical prediction for the finite size cor- 
rections of energies, the expected scaling was confirmed 
for both weak impurities and weak hopping. |^ However, 
due to the limited system size it was impossible to go be- 
yond the perturbative (in either Vb or 1^) regime. Later 
it was claimed that the full flow from a weak impurity 
to the open boundary flxed point (BFP) was successfully 
demonstrated, although this strong statement is 

not really supported by the numerical data presented. 
The smallest temperature discussed in Ref. [|2| corre- 
sponds to a system of around 300 lattice sites and the 
largest system considered in Ref. was A^ = 52, while 
in Ref. it was shown that A^ « 10^ lattice sites are 
clearly not enough to exclude an asymptotic behavior not 
governed by the BFP, even if one starts out with a fairly 
strong impurity. 

Recently functional RG methods, originally developed 
in a field theoretical context, have been introduced as a 
new powerful tool in the theory of interacting Fermi sys- 
tems [ p^ , with applications so far concentrating on trans- 
lationally invariant two-dimensional systems. p5[ | In this 
letter we apply such a functional RG scheme to the spin- 
less fermion model with site or hopping impurities. We 
compute the complete coupled flow of the renormalized 
onsite energies and the renormalized hopping amplitudes 
from the microscopic energy scale down to the infrared 
fixed point. The fiow equations are non-perturbative in 
the impurity strength while perturbative in the electron- 
electron interaction. We treat the full functional form of 
the renormalized impurity potential as generated by the 
flow, instead of replacing it approximately by the scatter- 
ing amplitudes at the Fermi level. Computing the local 
density of states near the impurity we convincingly con- 
firm the universality of the BFP. However, it turns out 
that very large systems (10^ — 10^ sites) are required to 
reach the BFP for intermediate impurity and interaction 
parameters. Our RG scheme is checked against numeri- 
cal exact DMRG data for systems with up to A^ = 768 
sites. 

The one-dimensional lattice model of spinless fermions 
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with nearest neighbor hopping ampHtude t = 1 and near- 
est neighbor interaction U is given by 



H = - 



+ Uy^njnj+i, (1) 



in standard second-quantized notation. Here we focus 
on the half filled band case. This model is either comple- 
mented by a site impurity Hs = Vrij^ or provided with 
a hopping impurity Hi^ — —tm{cj^c^^^^ + h.c.) on one of 
its bonds. 

In a weakly interacting spinless LL with an open end 
the local density of states Pj{uj) near the boundary can 
to a surprisingly good approximation be obtained from 
a non-selfconsistent Hartree-Fock (HF) approximation. 
Ip^ It is instructive to consider also the impurity problem 
within HF, before turning to the RG treatment. The im- 
purity leads to Friedel oscillations in the non-interacting 
density profile {nj)o which for large \j — jo| behaves as 
_Rsin {2kp\j — jo\)/\j — jo\, where R is the reflection am- 
plitude of the bare impurity. Similar oscillations are 
found in the matrix element {cjCj^-^)o. Thus both the 



only by a flnite amount of order U^. Hence, we can 
replace the renormalized two-particle vertex to leading 
order in U by the antisymmetrized bare interaction. In 
this way the exact hierarchy of flow equations gets trun- 
cated, and one obtains a simple one-loop flow equation 
for the self-energy E, where only the (full) propagator G 
and the bare electron-electron interaction U enter. Car- 
rying out a Matsubara sum and choosing a real space 
representation of G and E, one obtains the flow equa- 
tions (at temperature T = 0) 



dA 
dA^""±i 



—T^ = -— \^ \^ 

-JJ 27r ^j+sj+s 

s=±lw=±A 



(iuj) 



U 
2^ 



(3) 
(4) 



The self-energy is frequency independent and tridiagonal, 
since the bare interaction is instantaneous and restricted 
to nearest neighbors. The full propagator G^ on the right 
hand side of the flow equations is obtained by inverting 
the matrix [G°] — E^. The bare site and/or hopping 
Hartree potential J7((nj_i)o + (nj+i)o) and the Fock impurity enter as initial conditions for E'^ at A = oo. 



"hopping correction" are oscillating and of long range. 
One then has to solve a (non-trivial) one-particle prob- 
lem within such a potential and with modulated hopping. 
Taking into account the Hartree term only, the resulting 
spectral weight for — > shows power- law behavior 
with an exponent which is proportional to the amplitude 
UR of the oscillations. |0] We have checked numerically 
(for systems of up to 10° lattice sites p^ , p^ ) that this 
behavior is not changed when the Fock term is included. 
Thus due to the long range nature of the effective po- 
tential and the hopping modulation already HF yields a 
power-law for the spectral weight, but with an exponent 
which not only depends on U , but via R also on the bare 
impurity strength. 

It is tempting to extend the HF study using self- 
consistent HF. |19|] However, it turns out that an iter- 
ative solution of the selfconsistent HF equat ion le ads for 
all [/ to a charge density wave groundstate. 



For a site impurity V at jo, one sets '^f~j^ — V, and 
for a hopping impurity between jo and jo + 1, one has 



y^A— C30 y^A— cxj 

^ jo, 30 + 1 - ^ jo + 1,30 



-trn, while the other matrix ele- 



is qualitatively incorrect since a single impurity cannot 
change bulk properties of the system. 

We now treat the problem using a fermionic functional 
RG approach. Cutting off the free propagator on a scale 
A and differentiating with respect to this flow parame- 
ter, an exact infinite hierarchy of coupled differential fiow 
equations for the one-particle irreducible vertex functions 
can be derived. [p0|-p2t For the impurity problem it is 
technically advantageous to use a frequency cut-off for 
the free propagator 



G^'^iioj) = Q{\uj\ - K)G'^{iuj) 



where G'^ is the free propagator without cut-off and lo the 
Matsubara frequency. A fiows from oo to 0. For spinless 
fermions the electron-electron interaction is renormalized 



ments are initially zero. The above flow equations are 
non-perturbative in the impurity parameters, in contrast 
to the perturbative bosonic RG. Written in momentum 
space the different scattering channels E/j are coupled. 
The self-energy at A = can be given a simple phys- 
ical meaning: E^J° represents an effective one-particle 
potential and '^fj+i is an effective modulation of the 
hopping. To calculate Pj{i^) one determines the spectral 
weights of the remaining one-particle problem. 

For small impurity strength V, after transforming to 
momentum space and taking N — > oo, Eqs. (^ and (Q) 
can be solved analytically, as long as E'^ stays small. 
For the backscattering this gives E^^ -kp ^^^^^ 
7] = U[l — cos{2kp)]/{TrvF) and the Fermi velocity vp- 
19| which Xo leading order in U, the exponent 77 is just Kp — 1 
|]l6| which shows that the non-perturbative fermionic RG 
captures the power-law increase found in the perturbative 
bosonic RG. 

Numerically integrating the RG equations for finite 
systems we can go beyond the perturbative regime. In 
each step of the integration we have to invert an iV x 
matrix. If we assume open boundary conditions in Ho, 
[G*'] — E^ is tridiagonal in real space and the numeri- 
cal effort is considerably reduced. This allowed us to 
treat systems with up to 2^^ — 32768 lattice sites. For 
finite N the flow is effectively cut off on a scale of the 
order of For smaller systems we also considered 

periodic boundary conditions. Fig. |l| shows typical re- 
sults for E^J*' and S^J^j^ for a site impurity and lattice 
sites close to jo ■ Since E is symmetric around jo mainly 



(2) 
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the region j < jo is shown. Similar to HF the effective 
potential and hopping are oscillating and slowly decay- 
ing. The numerical data suggest that the oscillations 
again fall off as \j — jo\~^, but with an amplitude which 
compared to HF is slightly enhanced. The inset of Fig. 
|l| shows Sj^ as a function of A for different N. Ob- 
viously the renormalized potential at the impurity site 
remains finite and the expected "cutting" of the chain 
does certainly not occur because a single on-site energy 
diverges, as one might guess if the bosonic RG is taken 
too literally. Singular behavior is only found in j,, for 
momenta with fc — fc' w ±2kp, which is associated with 
the long range oscillations in real space. 
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FIG. 1 
jo = 513. 

(circles), A'' = 130 (squares), N = 258 (diamonds), A' 
(triangles), and N — 1026 (stars). 
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The inset shows 'Sj^jg as a function of A for A^ = 66 
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It shows E^^' , 1 as a 
We have checked that 



In Fig. g we present results for the case of a weak hop 
ping between two open chains vA 
function of A for = 1024. 
the curve to a good approximation already presents the 
A^ — > oo result. In contrast to a simplistic interpreta- 
tion of the bosonization result, the renormalized hopping 
io+i 'io^s '^''^ scale to zero. Similar to the case of a 
site impurity, S^J,° shows long range oscillations in both 
the effective potential and the hopping. Again this and 
not the scaling of a single Vj or tj is the reason for the 
peculiar behavior of physical observables, as for example 
Pj{Lo), discussed next. 

As an inset to Fig. 2 the density of states near the 
impurity, /3j,j_i(Lj), is presented for a site impurity. The 
data show a suppression of the weight for — > as 
expected. Each spike represents a (5-peak of the finite 
system. Instead of trying to fit a power-law to these 
data it is advantageous to analyze the finite size scaling 
of the spectral weight W{N) at ^. If pjo_i(w) fol- 
lows a power-law as a function of frequency, we expect a 
power-law with the same exponent in the A" dependence 
oiW{N). 



-0.10 



-0.11 3 



-0.12 




-0.13 -ooeeeeeoo^ 



U=1,ti=0.1 



10 ' 



10 



10 10 

A 



10 



10^ 



FIG. 2. Ejj,_j^_|_i as a function of A for a hopping impurity 
with tjg = 0.1, U = I, N = 1024, and jo = 512. Inset: 
Pjp_i(a;) as a function of uj for a site impurity with the same 
parameters as in Fig. 1 (A' = 1026). 



In Fig. H we show the negative of the logarithmic cen- 
tered differences ai{N) of W{N) as a function of A^ for 
U = 0.5 and different V obtained from RG (A^ < 32768) 
and DMRG {N < 768). If W{N) decays for A^ oo 
as a power-law, ai{N) converges to the respective expo- 
nent. For comparison we also calculated asiN) for the 
lattice site next to an open boundary (V — oo). The 
DMRG and RG data are parallel to each other, which in 
addition to the analytical arguments is a strong indica- 
tion that our fcrmionic RG captures the essential physics. 
For V — oo both methods produce the expected power- 
law behavior with boundary exponents a^'^^^ and a^*^ . 



,,DMRG 



(A^ — 512) agrees up to 1% with the exact expo- 
0.1609. a^'^iN = 16384) which effectively 



nent 

is equal to a^'^{N = oo) deviates by roughly 6% from 
ag^ since the RG is only correct to leading order in U. 
The RG curves for finite V suggest that for N ^ oo the 
af''^{N) converge to the universal {V independent) expo- 
nent a^*^ . This is in agreement with the field theoretical 
prediction. It is remarkable that even for fairly strong im- 
purities {V — 4) extremely large A^ = 10''-10^ are needed 
to exclude non-universal {V dependent) fixed points with 
some certainty. Solely relying on DMRG data for a few 
hundred lattice sites would in this case give no definite 
result. |l|] In Fig. | RG and DMRG data are presented 
for an intermediate impurity strength V = 1 and ^ = oo 
for different values of U. Due to higher order corrections 
in U, the difference between the RG and DMRG data 
increases with increasing U. For larger U, the af '^{N) 
approach faster, but even for the largest U = 1.8 
considered here ||2^ (which corresponds to Kp w 0.58) 
very large A^ are needed. This demonstrates that for 
intermediate V and U, which are experimentally most 
relevant, very large systems are needed to observe the 
universal BFP physics. For chains which are not long 
enough a strong system size dependence of experimen- 
tally extracted exponents must be expected. 
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FIG. 3. ai{N) as a function of iV for U = 0.5 and different 
V: V = 1 (circles), V = 2 (squares), V — 4 (diamonds), and 
V = oo (triangles). The filled symbols are DMRG data and 
the open ones obtained from the RG. The dashed-dotted line 
gives the exact boundary exponent a^. 
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FIG. 4. ai{N) as a function of TV for V = 1 (dashed 
lines) and as (N) for V = oo (solid lines) for different U : 
U = 0.5 (circles), U = 1 (squares), and U = 1.8 (diamonds). 
Filled symbols are DMRG data, open ones RG results. The 
dashed-dotted lines give the exact U dependent boundary ex- 
ponents a^g. 



We finally note that in the fermionic RG used in Ref. 
1^ flow equations were set up for a single parameter 
only: the transmission amplitude at the Fermi level. Our 
functional RG flow however indicates that in the non- 
perturbative regime different momentum channels are 
strongly coupled. Hence, we believe that it is impor- 
tant to take the whole renormalized impurity potential 
proflle into account. The RG equations used in Ref. ||] 
can also be derived within our formalism, if one makes 
similar crude approximations. [ p7| Also wc did not find 
signs of an enhanced spectral weight as predicted in Ref. 

In summary, by solving a functional flow equation in 
a fermionic representation we have shown that in a one- 
dimensional lattice electron system with Luttinger liquid 
behavior an impurity makes observables at low energy 



scales behave as if the chain is split in two parts with open 
boundary conditions at the end points. Our fermionic 
RG is non-perturbative in the impurity strength. Long- 
range oscillations in the effective impurity potential pro- 
vide a simple real-space picture of the "splitting" mech- 
anism. The accuracy of the finite site RG scheme was 
confirmed by a direct comparison to DMRG data. For re- 
alistic parameters very large systems are needed to reach 
the asymptotic open chain regime. Hence only special 
mesoscopic systems, such as very long carbon nanotubes, 
are suitable for experimentally observing the impurity- 
induced asymptotic open boundary physics. Our method 
can easily be generalized to the case of several impurities 
and e.g. resonance phenomena can be studied. ||l8|| 
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